Effect of topology on the critical charge in graphene 
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We show that the critical charge for the Dirac excitations in gapless graphene depends on the 
spatial topology of the sample. In particular, for graphene cones, the effective value of the critical 
charge can tend towards zero for a suitable angle of the conical sample. We discuss the nature of the 
scattering phase shifts, quasi-bound state energies and local density of states for a gapless graphene 
cone and determine the dependence of these physical quantities on the sample topology. 
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INTRODUCTION 

\ The fabrication of monolayer graphene [3-0| provides an opportunity to study the effect of sample topology on its 
electronic properties |4|- |9|. The low energy excitations of the gapless graphene are described by a two dimensional 
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' , ' massless Dirac equation [1CH14| . These excitations are negatively charged and behave like electrons with the Fermi 
velocity Vp ~ 10 6 m/s. In graphene, even a small external charge impurity Ze ~ 1 can produce strong nonperturbative 

,13 i effects, since the effective Coulomb interaction strength a = fc ~ 1 where the dielectric constant k ~ 5. Thus, 
graphene provides an ideal system to analyze the combined effects of topology and strong electric fields on its electronic 
D , properties. 

The idea of a critical charge plays an interesting role in this context U Tj^ . It is well known that massless Dirac 
fermions can tunnel through large external potentials, thereby avoiding trapping or formation of bound states. This 
phenomenon, known as Klein tunneling 2fj| - 24 1 , is possible in graphene due to the chiral nature of the Dirac excitations 
jl8ll20||. However, when the effective strength of the external charge exceeds a certain critical value a c , the low energy 
excitations can form quasi-bound states [15l - ll9j . For such states the wavefunctions oscillate rapidly near the position 
of the charge impurity, analogous to zitterbewegung in strong field QED 24, 25|. For planar graphene, the critical 
charge a c = ^ . The existence of such quasi-bound states and their effect on transport properties and local density of 
O . states (LDOS) in planar graphene have attracted considerable attention 15- ljj. 

In this paper we explore how a nontrivial topology of the sample may affect the critical charge and the corresponding 
, strong field QED effects in graphene. We address this issue for graphene cones, which are obtained by introducing 
' local defects in the hexagonal lattice structure, for example, in the form of a pentagonal ring |4|, [5J. The conical 
topology manifests itself in the form of nontrivial holonomies for the pseudoparticle wavefunctions. Such holonomies 
can be modelled through the introduction of a magnetic flux tube. Thus an electric charge localized at the apex of 
the cone can be equivalently described by a suitable combination of electric charge and magnetic flux tube. We show 
that the critical charge a c depends on the angle of the graphene cone. For certain angle of the graphene cone, the 
topological effects lead to the identification of the two Fermi points 0, HI ■ This leads to the surprising conclusion that 
for a certain value of the angle of the cone, the critical charge is zero. We discuss the effect of the conical topology 
T — I , on the scattering phase shifts, quasi-bound state energies and local density of states (LDOS) in graphene. 

This paper is organized as follows. In the next Section we set up the Dirac equation for gapless graphene cone 
. , with a point charge at the apex. This is followed by the analysis of the corresponding spectrum, where we obtain the 
scattering phase shits, quasi-bound state energies and local density of states (LDOS) and show how these physical 



quantities depend explicitly on the sample topology. We also briefly comment on the RG flow of the strength of the 
charge impurity. We end this paper with some discussion and outlook. 
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DIRAC EQUATION FOR A GRAPHENE CONE WITH A COULOMB CHARGE 

In this section we set up the Dirac equation for a graphene cone with a coulomb charge at the apex of the cone. 
The low energy properties of the quasip article states near the Dirac points in graphene can be described by the four 
component Dirac wave function lB- 12j given by 



The pseudospin indices A and B label the two sublattices of the primitive cell of graphene and the valley indices + 
and — label the two inequivalent Dirac points K+ and K_ respectively. 

Considering low energy excitations about the Dirac point K + , Dirac equation for a planar gapless graphene in the 
presence of a Coulomb charge Ze is given by 



(X 
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# = EV, (2) 



where r is the radial coordinate in the two dimensional x — y plane. The Pauli matrices <Ji,2,3 and the identity matrix 
<7o act on the pseudospin indices^, B. The Coulomb interaction strength a = jf K e Vp ■ From now on we set H = vp = 1. 
In a two dimensional plane, the angular boundary condition for a Dirac spinor as it goes around a closed path is given 
by 

tf(r,0 = 2tt) = e 47r<T3 *(r,6» = 0). (3) 

A cone is obtained from the two dimensional plane by introducing a topological defect, which modifies the angular 
boundary condition ([3]). Following [J], consider a cone, formed by removing the sector AOB (see Fig. la.) and then 
identifying the edge OA, labeled by 8 = with the edge OB, labeled by 8 = 2ir. The sector AOB subtends an angle 
^ at the centre O. Due to this identification the frame {e x , e y } becomes discontinuous across the joining line. This 
problem can be solved by choosing a new set of frames given by 

e x > = e e and ey = -e r , (4) 

which is rotated with respect to old frame (e x , e y ) by an angle 4> = 8+^ in the counter clockwise direction, see Fig. lb. 
To keep the form of the Hamiltonian the same, the wave function has to be transformed by exp(i</>0"3/2). Similarly, if 
n sectors are removed from the plane where n can take only discrete values 1, 2, 3, 4, 5, the wave function should be 
transformed as 

*(r, 8 = 2n) = -e i2n(1 -^^{r, 8 = 0). (5) 

Note that an additional negative sign appears in the boundary condition since we take the angular part of the 
wavefunction to be e^ e , j being a half integer 0, [B| . 

When n is odd, an additional phase change is required due to the breaking of the bipartite nature of the hexagonal 
lattice. It can be seen from the Fig.((2|), that the removal of a wedge of opening angle introduces a pentagonal 
defect. When the cone is formed by removing a single wedge of angle the adjacent points on two sides of the 
identification line belong to the same sublattice. This feature is valid for all allowed odd n = 1,3,5 and leads to a 
system with only one Fermi point. Considering all the factors, the angular boundary condition satisfied for a general 
value of n is given by [J, [5j . 

*(r,6» = 27r) = -e i27r I- Ii ?+( 1 -f^]*(r,6' = 0). (6) 

Here Ti acts on the valley indices (+, — ). Note that for n = 0, the angular boundary condition for the planar graphene 
is recovered. For the allowed even values of n = 2 and 4, it can be seen from © that the matrix t 2 plays no role in the 
angular boundary condition This is due to the fact that for even n, e 12 ^^] = cos(^y 5: ) = e Z7r I~l. For odd 

n, the two Dirac points are identified and the energy eigenstates of the Dirac equation are obtained by diagonalizing 
T2 which mixes the valley indices in the wavefunction [4, 5]. Since the eigenvalues of T2 are ±1, for all values of n, we 
can replace the angular boundary condition ([6]) by 

*(r, 8 = 2tt) = _ e ^[±f +(i-f Wyfag = ). (7) 
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FIG. 1: (a)Formation of a cone from plane graphene sheet by cut and paste procedure and (b)Rotation of the coordinate frame 
due to its new orientation. In Fig. (a) blue atoms represent sublattice A and red atoms represent sublattice B. 




FIG. 2: A Graphene cone formed by removing an angular sector of -jr from the plane sheet. 



In subsequent analysis, we shall use the form ([7} of the angular boundary condition assuming that for all allowed odd 
values of n, the matrix T2 acting on the valley indices has already been diagonalized. 

The effect of the angular boundary condition ([7} on the wave function can be equivalently described by introducing 
a magnetic flux tube passing through the apex of the cone The presence of a magnetic vector potential modifies 
the boundary condition on a Dirac spinor as 



*(r,0 = 2tt) = -e te * Xa *(r,0 = 0). 

Taking dl as a line element on the circumference of the cone at a distance r from the apex, we have 

-* n 
dl = e g r(l - -)d6. 

Substituting © in © and assuming that A s is independent of the angle 9, we get 

*(r, = 2tt) = _ e te2 '"-( 1 -S ) A «*(r, = 0). 



(8) 



(9) 



(10) 
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TABLE I: The values of critical charge ot c i.e the minimum values of \u\ for different values of opening angle of the graphene 
cone i.e for different values of n. The values of angular momentum j for which we get the critical charge are also noted. 



value of n 


Critical charge (a c ) 


Corresponding angular momentum (j) 





0.5 




1 


0.3 


±1 


2 





±1 


3 


0.5 


±1 


4 


1.5 


±k 


5 


1.5 


±1 



Comparing equation ([7J) and equation (|10jl we get 

1 — 

Thus the entire effect of the conical topology on the wavefunction can be described by introducing the vector 
potential (10), which replaces the ordinary derivatives by the corresponding covariant derivatives in the Hamiltonian 
(|2|) . In this effective description, the Hamiltonian for a graphene cone with a Coulomb charge at the apex is given by 

H = -i{oid x > + o 2 dy>) - e(a ■ e g )A e - a (jy^j > ( 12 ) 

where the relation between (x',y') and (r, 9) can be obtained from the Equation^, which gives d x i — j._ ai d$ and 
d y i = —d r and a ■ eg = <J\. Thus the Dirac equation ([2]) can be written as 

We use an ansatz for the wavefunction given by 

where the total angular momentum j takes all half integer values. Substituting (fT4"| in dT51) . we note that the leading 
short distance behaviour of the wavefunction is given by 

*Jf B (r) ~ r 7 "2 where 7 = v 7 ^ 2 - a 2 and i/ = T^fv - (15) 

(1 — G' 

From (I15p it follows that when \a\ exceeds \u\, 7 becomes imaginary. As a result, the eigenstates ^>^\r) and 
(r) becomes wildly oscillatory and have no well defined limit as r — > 0, which corresponds to the phenomenon of 
zitterbewegung in a strong Coulomb electric field 25] . The critical value of the coupling is denoted by a c and it is 
given by the minimum allowed value of \v\. The parameter v depends on j and the number of sectors n removed from 
a plane to form the graphene cone. Hence we see that the critical coupling a c explicitly depends on the angle of the 
graphene cone. The values of the critical charge for different values of n are given in TABLE [J] 

From these values of critical charge, a surprising result can be observed for the case of n = 2, where a c is zero. In 
this case, any external charge will be supercritical. We should also note that when n = 5, none of the values of v in 
the lowest angular momentum channel j = corresponds to the minimum value of \u\ or the critical charge. In 
this case the critical charge of the system is obtained from the angular momentum channel j — ±|. 



SPECTRUM IN GRAPHENE CONE WITH A SUPERCRITICAL CHARGE 



In this Section we derive the scattering and bound state spectra for a gapless graphene cone in the presence of a 
supercritical Coulomb charge located at the apex of the cone. First consider the case v 7^ 0. Using (fT5|) and (|T4"1) . the 
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Dirac equation in each angular momentum channel j can be written as 



{dr + (v + m\ ( *<?(r) 



i*g>(r) 



111 terms of the functions u^\r) and w^'(r) defined by (r) = [w^'(r) + u^'^(r)] and (r) = [t)W)(r) 



e r ' 

1 r, n A VTV 0), 



(16) 



we get 



and 



r dV - — + (7 + (r) - i/u^ (r) = 



( r ) ... ... 

r r-^ + (7 - ia - 2iEr)u^> (r) - i/t)W' (r) = 0. 



From (H7J) and fig]), we get 



s r^r-^ + (1 + 27- s) — - (7 + ia)f U) (s) = 0, 



ds 2 



r/.s 



(17) 



(18) 



(19) 



where s — ~2ikr, with k = —E. For the discussion below, for any given value of n and j, we choose a greater than 
the corresponding value of \v\. This ensures that the coupling is always in the supercritical region. Define 7 = iX 
where A = y/a 2 — v 2 . Then the solution of Eq. (|19p is given by 

( s ) = dM (i(A + a), 1 + 2£A, s) + C 2 s~ 2a M (z(a - A), 1 - 2iA, s) . (20) 

From ((TTJ) and ([2D]) we get 

u (i) (s) = -idfiM (1 + i(A + a), 1 + 2iA, s) - i(C 2 / ix)s- 2iX M (1 + i(a — A), 1 - 2iA, s) (21) 



where fi = yg^j. 

In order to proceed, we use the zigzag edge boundary condition [u^\ao) — vV'(oq)] = 0, where cto is a distance 
from the apex, of the order of the lattice scale in graphene [16| ■ This gives 



C 2 = e 2l ^ k) ue* x Ci where 



,2c 



C(fc) _ »(1 + ^) 2iAln(2fcq ) 



(1 - i M ) 

From the above, we obtain the scattering matrix S as 



£ _ e 2i5„(fe) 



2ialn(2fcr) 



where f a ,\= rjl+iA-^i) • F rom Q23[) we obtain the scattering phase as 



$„(*:) = arg[e- 4C(fc) + be lC{k) ] - a\n{2kr) + arg(/ QjA ) 



(22) 



(23) 



(24) 



where b — e 7rA /j ^J' ^ A . In Fig. (3a) we plot (fM|) ignoring the Coulomb tail term —aln(2kr). For ^ = 0.3, the value of 
a = 1.8 is deeper in the supercritical region than it is for v = 1.5. When the coupling a is deeper in the supercritical 
region, the phase shift has more number of kinks, which indicate the bound states. The plot also shows that the phase 
shift depends on the topology through its dependence on n via v. 

As mentioned before, in gapless graphene we do not expect bound states due to Klein tunneling. However, in 
the supercritical regime, the system admits quasi-bound states whose spectrum is obtained from the zeroes of the S 
matrix in (|23[) . The quasi-bound state energies are given by 



kp = i eXP 



1 

2^ 



(25) 



where p is a positive integer. 
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FIG. 3: (3a) shows the dependence of scattering phase shift 6 on wavenumber kao for v — 0.3, 0.5, 0.9, 1.5 and a = 1.8, ignoring 
the Coulomb tail term —aln(2kr). As the value of v increases, the kinks in the phase shift become sharper, which indicates 
the dependence of the phase shift on the angle of the graphene cone. (3b) shows dependence of ground state energy on the 

3 + — "I 

Coulomb potential strength for different angles of the graphene cone. We have considered v — 1 _ t and j = 5. 




Wavenumber (kr) Wavenumber (kr) 

FIG. 4: Energy dependence of LDOS in presence of a Coulomb potential for different values of n and a particular value of r 
and with j = i and a = 0.6. This is in the supercritical region for both n — and 1. 



Another interesting observable in this context is the LDOS. In Fig.(Q| we have plotted the standing wave oscillations 
in LDOS p(k, r) using 

P(fc,r) = -i-^|*^(fc,r)| 2 , (26) 

3 

where (k, r) is the radial part of the spinor given in Eq. (|14[) . for a given angular momentum channel j. The maxima 
and minima of the oscillations occur at half-integral and integral multiple of n respectively. The normalization constant 
of the wave function depends on the value of n which decreases with increasing n. Therefore, though the nature of 
standing wave oscillation is similar for n = and 1, the amplitudes of the corresponding LDOS decreases with the 
increase of n. 

It is known that in the planar case, under the RG flow, a supercritical charge is driven to its critical value [l6j |. 
While the full discussion of the RG flow requires a more detailed analysis, we can obtain some of its qualitative 
features from the quasi-bound state energies. To see this consider the real part of k p in (|25|) which gives the energy. 
This part diverges as the cutoff ap — > 0. In order to study the RG flow, we now promote the coupling constant a as 
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a function of ao and demand that as ao — > 0, the energy for any fixed level p (say p = 1) remains independent of the 
cutoff 2(j], 27 1 - In the leading order, where a is only slightly above the critical coupling, this prescription gives the 



/3-function as 

/3(A)- -A 2 + .. (27) 

Thus we have an ultraviolet stable fixed point at A = or at a = v. Hence, for any given value of n and j, the 
coupling a in the supercritical regime is driven to its critical value. 

We end this Section with a comment for the case v = 0, which occurs for n = 2 and j = ±i. For v — 0, the 
first order Dirac equations (|17[) and (I18p completely decouple and the functions u and v become independent of each 
other. In this case, the critical coupling a c also vanishes. This means that any external charge, no matter how small, 
would lead to the supercritical regime. 



CONCLUSION 



In this paper we use the fact that system of a graphene cone with an external Coulomb charge at its apex can be 
equivalently described by the combination of the Coulomb charge and a suitable magnetic flux tube passing through 
the apex. The above analysis has been done in the supercritical regime, where the external Coulomb charge exceeds 
a certain critical value. 

The quantities of physical interest in this system include the scattering phase shifts, the LDOS and the quasi-bound 
state energies. We have shown that all these physical quantities depend explicitly on the number of sectors removed 
from a planar graphene to form the cone. The existence of the quasi-bound states indicates the possibility of the 
localization of the wavefunctions in the presence of a supercritical charge. Our analysis shows that the nature and 
extent of the localization depends on the spatial topology of the graphene sample. 

Finally, we have given qualitative arguments which shows that under the RG flow and for v =/= 0, the supercritical 
charge in the graphene cone tends to its critical value. If this argument can be extended for v = 0, for which the 
critical charge vanishes, that would lead to complete shielding of the external charge. This issue and the related 
electronic properties are currently under investigation. 
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